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ABSTRACT 

The velocity anisotropy of particles inside dark matter (DM) haloes is an important 
physical quantity, which is required for the accurate modelling of mass profiles of 
galaxies and clusters of galaxies. It is typically measured using the ratio of the radial- 
to-tangential velocity dispersions at a given distance from the halo centre. However, 
this measure is insufficient to describe the dynamics of realistic haloes, which are 
not spherical and are typically quite elongated. Studying the velocity distribution 
in massive DM haloes in cosmological simulations, we find that in the inner parts 
of the haloes the local velocity ellipsoids are strongly aligned with the major axis 
of the halo, the alignment being stronger for more relaxed haloes. In the outer re- 
gions of the haloes, the alignment becomes gradually weaker and the orientation is 
more random. These two distinct regions of different degree of the alignment coincide 
with two characteristic regimes of the DM density profile: a shallow inner cusp and a 
steep outer profile that are separated by a characteristic radius at which the density 
declines as p oc r -2 . This alignment of the local velocity ellipsoids requires reinterpre- 
tation of features found in measurements based on the spherically averaged ratio of 
the radial-to-tangential velocity dispersions. In particular, we show that the velocity 
distribution in the central halo regions is highly anisotropic. For cluster-size haloes 
with mass 10 14 — 10 15 /i _1 M Q , the velocity anisotropy along the major axis is nearly 
independent of radius and is equal to (3 = 1 — Cperp/^radial ~ 0.4, which is significantly 
larger than the previously estimated spherically averaged velocity anisotropy. The 
alignment of density and velocity anisotropics, and the radial trends may also have 
some implications for the mass modelling based on kinematical data of such objects 
as galaxy clusters or dwarf spheroidals, where the orbital anisotropy is a key element 
in an unbiased mass inference. 

Key words: galaxies: clusters: general - galaxies: kinematics and dynamics - cos- 
mology: dark matter 



1 INTRODUCTION 

The orbital anisotropy describes the distribution of or- 
bits in astrophysical systems. It plays a key role 
in dynamical modelling of kinematical data o f ob- 
jects at all scales , from dwarf spheroi dals dLokasI 



Walker ct al. 2009), to elliptical galaxies (iDekel et alj 



2009 



2005 



Napolitano et al.ll201l1:IWoitak fc Mamonll2012l) ancTgalax y 
clusters (|Biviano fc Girardil 120031 : IWoitak fc Lokasl l2010l) . 
Prior knowledge on the anisotropy or elaborated techniques 
of data analysis are essential for accurate and unbiased 
mass estimates. The main difficulty arises from the well- 
known mass-anisotropy degeneracy occurri ng in the Jean s 
analysis of the velocity dispersion profiles (|Merrittl 1 19871 ). 
Several methods were developed to break this degeneracy 



(see e.g. lLokasi |2002|; lLokas fc Mamonl 120031 : IWoitak et all 
l2009l : IWolf et al.ll2'oTol : lMamon et al.ll2012f TTheir efficiency, 
however, critically relies on the quality of the data. In ad- 
dition, observational constraints on the mass profiles are 
still affected by what dynamical models assume about the 
anisotropy or what priors on the anisotropy are used. As 
the matter of fact, most mo dels rely on certain parametrisa- 
tion of the anisotro py (e.g. Hookas 2002; M amon et al.ll2012l : 
IWoitak et al.ll200^) or incorporate some we ll-motivated pro- 
files (e.g. iDiaferid [l999l : IDekel et al.1l2005h or prior distri- 
butions for its parameters (e.g. iNewman et al.l l2012l ). The 
choice of the prior probability and parametrisation is com- 
monly motivated by cosmological iV-body simulations. The 
orbital anisotropy in DM haloes is often used as a point of 
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reference in such preselection of dynamical models. The best 
example showing this effect of feedback from the simulations 
is commonly adopted assumption that the velocity distri- 
bution in the centre of gravitationally bound cosmological 
objects is isotropic or near ly isotropic (see e.g. IWalker et al.l 
I2006I : In cwman et al.ll20ll as examples in studies of dwarf 
spheroidals and galaxy clusters). 

The orbital anisotropy is also one of the key parameters 
describing the phase-space structure of DM haloes. Many 
studies showed a number of interesting properties such as 
a relation between the anisotropy and D M density profile 
( Hansen fc Moord liool IZait et al.ll2008l l or existence of a 
universal attractor in t he space spanned by all solutions 
of the Jeans equation ( Hansen et al.l |2010f) . As shown by 
iDehnen fc McLaughlin] (|2005l ). the anisotropy and its rela- 
tion to the slope of the density profile is also a key element in 
theoretical understanding of the universal phase-space den- 
sity profiles measured in DM haloes. A number of studies ad- 
dressed such problems as the radial profiles of the a nisotropy 
(|Woitak et al.ll2005l ; lAscasibar fc Gottloberll200cf), the bia s 
between DM particles and subhaloes (|Diemand et al.ll2004h . 
evolution of the anisotropy in con trolled simulations of 
halo mergers dSparre fc Hanserj|2012l ), the redshift evolution 
ijlannuzzi fc Dolaell2012| ) and dependa nce on the mass an d 
dynamical equilibrium of DM haloes (|Lemze et al.l I20I2I ) . 
The overall picture emerging from these studies points to 
the fact that the velocity distribution is nearly isotropic in 
the halo centre and radially biased at large radii (excess of 
radial orbits). Although this trend is common to all DM 
haloes, the anisotropy profiles of individual haloes are sig- 
nificantly scattered around the mean trend. Deviation from 
the me an trend seems to depend on dynamical stage of the 
haloes (jLemze et al.l[2012l ). 

The global velocity ellipsoids of DM haloes are aligned 
with the major axes of the halo shape (Kasun fc Evrardl 
2003; lAlkood et al 1 120061 : ISaro et al.ll2012r T This property 
seems to occur also locally, although thi s problem has no t 
been extensively studied in the literature (|Zemp et al .120091) . 
Interestingly, similar configuration of the velocity ellipsoid 
is consistent with the k inematical data of ell iptical galaxies 
llCappellari et al.ll2007l ) and galaxy clusters (|Skielboe et alj 
l2012t ) what suggests that it is probably a generic feature 
of the orbital structure in cosmological objects. Despite 
these facts, the orbital anisotropy in simulated DM haloes 
is commonly quantified in terms of the ratio of the radial- 
to-tangential velocity dispersions which breaks the preferred 
symmetry of the velocity distributions. This raises the ques- 
tion to what extent the orbital anisotropy defined in this 
way describes the true distribution of orbits in DM haloes. 

The anisotropy of the velocity dispersion tensor 
is commonly quantified by the anisotropy parameter 
l|Binnev fc Tremaindl200sl ) 



P 2ol ' 



(1) 



where of and of are the radial and tangential velocity dis- 
persions, respectively. This is a reasonable definition for a 
spherical or almost spherical halo because in this case the 
velocity dispersion tensor should be oriented along the ra- 
dial direction and, because of the symmetry, one expects 
that dispersions perpendicular to the radius are equal: a'g — 
a\ = cf/2. For haloes with substantial elongation - and a 



large fraction of cluster-size haloes are elongated - the veloc- 
ity ellipsoid is more complex and simple spherically averaged 
/3(r) is not enough to describe such systems. 

Despite obvious inconsistencies between the spherically 
averaged j3 and the elongation of DM haloes, the anisotropy 
parameter f3 is commonly used to study the phase-space 
properties of DM haloes. In realistic cases one expects that 
local velocity tensor is neither aligned along radius nor tan- 
gential velocity dispersions are equal. To complicate the sit- 
uation, the velocity dispersions are also expected to depend 
on all three space coordinates, not just one radius. So, real- 
istic /3(r) is a complex entity. To some degree, this justifies 
using spherically averaged /3. Realistic /3(f) is not well stud- 
ied and, even it were, it would be more difficult to adopt it 
in dynamical models. Instead, by using spherically averaged 
ft, one gets a practical tool to build dynamical models with 
understanding that there is a price for this simplification: 
one should expect some errors in the models. Here we do 
not study the errors and focus only on description of the 
realistic velocity ellipsoids. 

In this paper, we study the properties of the velocity 
ellipsoids in cluster-size DM haloes in high-resolution cos- 
mological simulations. The paper is organised as follows. 
In section 2, we describe cosmological simulations and the 
properties of DM haloes used for the analysis. In section 3, 
we investigate the local values of the anisotropy parameter 
defined using the radial and tangential velocity dispersions. 
This gives us insights into dynamical structure of the halo 
and it shows that anisotropy parameter depends on the po- 
sition with respect to the halo major axis. This "classical" 
treatment of the anisotropy parameter neglects the fact that 
velocity ellipsoid may not be aligned along radius and that 
it may be triaxial. In section 4, we study the properties of 
the local velocity dispersion tensor, e.g. alignment, triaxi- 
ality and anisotropy, and we show that the local velocity 
ellipsoids are aligned with the major axis of the halo shape. 
In section 5, we calculate spherically averaged (measured 
in spherical shells) profiles of the orbital anisotropy defined 
in the symmetry consistent with the alignment of the local 
velocity ellipsoids and compare with the local values. We 
also show how the anisotropy parameter based on the ratio 
of the radial-to-tangential velocity dispersions misrepresents 
the true orbital structure in DM haloes. Section 6 contains 
discussion and conclusions. In this section, we also provide 
a simple analytical model for the local velocity ellipsoids 
measured in DM haloes. 



2 METHOD 
2.1 Simulations 

We use cluste r-size DM haloes s elected from the Bolshoi 
simulation Q ilKlvpin et alj l201ll ) . The simulation follows 
the evolution of DM structures in the framework defined 
by a ACDM cosmological model with cosmological parame- 
ters consistent with r ecent measurements b ased on WMAP 
five-year data release (jKomatsu et al.ll2009l ) and abundance 



The simulation is publicly available t hrough the MultiD ark 
database ( |http: / /www.mul tidark.org). See lRiebe et al.N201ll) for 
all details of the database. 



© 0000 RAS, MNRAS 000, 000-000 



Orbital anisotropy 3 



r=0.05r. 




* * 



Figure 1. Density and velocity anisotropy of DM haloes as a function of the polar and azimuthal angles 8 and <f>. The system of 
coordinates is chosen in such away that the major axis is placed at 8 = tt/2 and <j> = 7r±7r/2. Left panels: Deviations of the local density 
from the average density at given radius p(r, 8, <f>) / p(r). Right panels: Local velocity anisotropy parameter 0. The rows show the maps 
calculated inside three spherical shells of radii equal to 0.05, 0.20 and 1.00 of the virial radius. When measured in spherical shells the 
density anisotropy declines with decreasing radius. However, the pattern of elongation along the major axis is clearly preserved at all 
radii. The velocity anisotropy 0(r,9,(j>) behaves very differently as compared with the density p(r,8,<p). Close to the virial radius the 
velocities are predominantly radial in all directions: m 0.4. At smaller distances has same value ~ 0.4 along the major axis and gets 
to < in the plain perpendicular to the major axis. 



of the Sloan Digital Sky Survey clusters ijRozo et al.ll201oT ). 
The simulation box has the size of 250/i _1 Mpc and con- 
tains 2048 3 particles, each with mass of 1.35 x 10 8 /i _1 M o . 
High mass resolution of the simulation is a key feature to 
study the properties of the local velocity distributions in DM 
haloes. Selected DM haloes contain from 5 x 10 5 to 8 x 10 6 
particles inside the virial spher e. For all deta i ls on the sim- 
ulations, we refer the reader to iKlypin et all l|201lf ). 



DM haloes are found using the Bound-De nsity-Maxima 
(BDM) algorithm jKlypin & Holtzmanl [19971 ). We use only 
distinct haloes and do not include subhaloes. BDM finds 
all density maxima with density estimated with the top- 
hat filter containing 20 particles. Among all density maxima 
inside a given distinct halo the code finds the one, which 
has the deepest gravitational potential, and uses it as the 
centre of the halo. The halo virial mass M v is defined as 
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a spherical overdensity mass with the mean density A = 
3M v /(47rrJ) = 97.2 times greater than the critical density, 
where r v is the virial radius. For our study, we select all 517 
haloes with virial masses greater than 10 14 Mq. Velocities 
of DM particles are corrected for the halo bulk velocities 
approximated by the mean velocity of the particles inside 
the virial sphere. 

2.2 Halo shapes 

We quantify the shape of the haloes in terms of the moment 
of inertia tensor given by 

N 

— ^ ] Tj,nTj,n, (2) 
n=l 

where n t n is i-th component of the position vector with re- 
spect to the halo centre and the sum is over all particles 
lying inside the virial sphere. Eigenvectors of the tensor de- 
termine the principle axes of an ellipsoid approximating the 
halo shape. Information which is essential for our study is 
the eigenvector associated with the major axis of the halo 
shape ellipsoid (the axis minimising the moment of inertia). 

2.3 Relaxed haloes 

In our analysis, we also consider a subsample of relaxed 
haloes. Our criteria of rel axness are similar to those pro- 
posed bv lNeto et all (|2007f ) and are based on three diagnos- 
tics of dynamical equilibrium: the offset between the mass 
centre and the minimum of the potential, the offset between 
the bulk velocity of the halo and the mean velocity of the 
most gravitationally bound particles, and the virial ratio. 
We use BDM halo centres as the positions of the minimum 
of the potential and BDM bulk velocities as velocities of the 
most gravitationally bound particles. The virial ratio is es- 
timated using 3 x 10 J particles randomly drawn from every 
halo. The gravitational binding energy is computed using 
direct summation. 

We define relaxed haloes as those which satisfy three 
conditions: the offset of the mass centre less than 0.07r„, the 
offset of the bulk velocity less than 0.2V V and the virial ratio 
less than 1.35. The limits imposed on all diagnostics separate 
the outliers which populate long tails of the distributions and 
are associated with unrelaxed haloes. Most unrelaxed haloes 
(76 per cent) are found by the offset of the mass centre. Our 
selection criteria yield 267 relaxed haloes (52 per cent of the 
total number). 



3 VELOCITY ANISOTROPY ALONG AND 
PERPENDICULAR TO THE MAJOR AXIS 

We start our analysis by investigating the global trends in 
the velocity anisotropy. Here we ask a simple question: how 
does P change along and perpendicular to the major axis of 
the density distribution. 

We compute the local values of the anisotropy param- 
eter on a regular grid of spherical angles inside spherical 
shells. The polar 6 and azimuthal cj> angles of the halo ma- 
jor axis are set consistently at the same values in all haloes 
(for the sake of readability of the plots, we chose 6 = n/2 
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Figure 2. Radial profiles of the anisotropy parameter /3. The pro- 
files are complementary to Figure [l] showing dependance of /3 on 
the position with respect to the halo major axis. Spherically av- 
eraged /3(r) indicates preferentially radial anisotropy at all radii. 
It also slightly increases with radius, as found in many studies. 
Instead, the anisotropy along the density major axis is nearly con- 
stant at all radii (with /3 Rj 0.3), whereas regions perpendicular to 
the major axis show tendency to have almost isotropic velocities 
at small distances (/? < 0). Shaded regions on the plot show 50 
per cent statistical uncertainties. 

and (j> — 7r ± 7r/2). We use three spherical shells of radii 
(0.03 - 0.07)r v , (0.16 - 0.25)r v and (0.93 - 1.07)r v contain- 
ing on average 3.9 x 10 4 , 1.2 x 10 5 and 10 5 particles. The 
two inner shells probe the inner region of the halo where 
the spherically averaged DM density profile changes grad- 
ually from p(r) cx r _1 to p(r) cx r~ 2 , and the third one 
is the virial sphere. The anisotropy parameter is calculated 
inside the cones with opening angle of 10 degrees around ev- 
ery point of the regular grid representing a cylindrical map 
projection of a sphere (the Mercator projection), for every 
spherical shell. At every grid point, we find the median value 
of the anisotropy parameter in the halo sample. In all results 
presented below we show the medial values. 

The right panels of Figure [T] show the dependance of 
the anisotropy parameter /3 on the position with respect to 
the halo major axis. The left panels show the values of the 
ratio p(r, 9, <j>)/ p(r), where p(r, 8, (j>) is the density inside the 
intersection of a radial shell and a cone directed to (6,cj>), 
and p(r) is the density inside a shell. 

The maps of the local anisotropy parameter inside the 
two inner shells show a dipole-like dependance on the po- 
sition with respect to the halo major axis. The value of /3 
changes from « 0.3 along the major axis to fs in the plane 
perpendicular to the major axis. Similar dipole structure, 
although not as prominent as in the inner shells, is also visi- 
ble at the virial sphere. Here, the local anisotropy parameter 
tends to change from 0.5 (along major axis) to 0.3 (in the 
perpendicular plane). We note that the same spatial vari- 
ation of the local anisotropy parameter has been recently 
found in non-c osmological controlled si mulations of mergers 
of DM haloes (jSparre fc Hansen! 120121 '). In every shell, the 
dipole of the anisotropy parameter coincides with the dipole 
of the density ratio (compare the left and right panels in 
Figure [TJ . 

The values of the spherically averaged velocity 
anisotropy (computed in spherical shells) are 0.14, 0.20 and 
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0.30, from the innermost to outermost shell respectively. As 
Figure [2] shows, spherically averaged values of /3 are substan- 
tially different from /3 along and perpendicular to the major 
halo axis. This discrepancy between the local and spheri- 
cally averaged parameters is the most prominent in the two 
inner shells. 

Alignment of the velocity anisotropy dipole with the 
density dipole suggests that the apparent dependance of /3 
on the position with respect to the halo major axis may be a 
purely geometrical effect resulting from a wrong assumption 
that the local velocity ellipsoids are aligned with the radial 
directions (assumption underlying the definition of /3). In 
order to address this problem in more detail, we shall study 
the properties of the local velocity dispersion tensors. 



4 LOCAL VELOCITY ELLIPSOIDS 

The anisotropy parameter ,3 = 1 — o 2 /(2of ) measures the 
true anisotropy of the velocity distribution only in spheri- 
cal systems, where one can assume that velocity dispersion 
tensor is diagonal in spherical coordinate system and two 
diagonal coefficients are equal: of — of = of/2. In non- 
sphcrical systems with distinct centres such as DM haloes, /3 
is just a function of 6 non-vanishing coefficients of the veloc- 
ity dispersion tensor. This makes this parameter incapable 
of disentangling the anisotropy of the velocity distribution 
from its spatial orientation (shape from orientation of the 
local velocity ellipsoids - geometrical representations of the 
velocity dispersion tensor). In particular, ft = may result 
from non-radially oriented and elongated velocity ellipsoids. 

In this section, we measure all components of the ve- 
locity dispersion tensor at different positions in DM haloes. 
Complete information on the tensor allows to disentangle 
the shape from the orientation of the local velocity ellip- 
soids. The local velocity dispersion tensor is calculated in a 
volume element around a certain position as 

1 N 

= — ^ y ] Vi,nVj : n (3) 
n — 1 

where Vi are the velocity components of the particles with 
respect to the mean velocity of the volume element and N 
is the number of DM particles inside the volume element. 

We compute the velocity dispersion tensor in several 
spherical shells on a grid deter mined by sphere p ixelisation 
defined by the HEALPix code (|G6rski et alj|2005h . with the 
total number of pixels fixed at 48. The polar angle of the 
pixels is measured with respect to the halo major axis. This 
allows for an easy selection of the pixels at equal angular 
distances from the halo major axis. We consider two sub- 
samples of the pixels: 8 pixels around the major axis with 
cos#| > 2/3 and 8 pixels in the plane perpendicular to the 
major axis, with |cos#| < 1/3. 

We select DM particles using six spherical shells of 
radii: (0.008 - 0.015)r V) (0.04 - 0.06)r v , (0.09 - 0.11)r v , 
(0.19 - 0.21)rv, (0.48 - 0.52)r v and (0.98 - 1.02)r v . This 
choice of the shells leads to 10 2 — 10 3 particles per pixel 
(and 20 — 100 for the innermost shell). All six independent 
components of the velocity dispersion tensor are calculated 
using the eq. Q. We diagonalise the local velocity disper- 
sion tensor obtaining three eigenvalues of > of > of and 
three associated eigenvectors. 



We calculate the distributions of various properties of 
the local velocity dispersion tensor in the halo sample. The 
distributions are computed by combining information from 
the same sets of pixels in all DM haloes. This means that the 
resulting distributions represent pixel-weighted statistics (all 
haloes contribute equally to the distributions). We consider 
several combinations of the pixel and halo subsamples: all 
48 pixels from all haloes (48 x 517 data points), 8 pixels 
along the major axis (or in the perpendicular plane) from 
all haloes (8 x 517 data points), 48 pixels from relaxed haloes 
(48 x 267 data points). 

4.1 Prolateness 

Geometrical representation of the velocity dispersion tensor 
is a triaxial ellipsoid with the axes proportional to oi, 02 
and 03. The shape of such ellipsoid may be characterised by 
the triaxiality parameter 

Z = (4) 
02 — 03 

introduced by iBinnevI (|l985l ) for description of the shape 
in the position space and adopted here for the veloc- 
ity space. This parameter measures the degree of prolate- 
ness/oblateness of the ellipsoids. Two limiting cases with 
Z = and Z = 1 correspond to oblate and prolate ellip- 
soids, respectively. 

Figure [3] shows the distributions of the triaxiality pa- 
rameter Z of the local velocity ellipsoids. The ellipsoids tend 
to be preferentially prolate at all radii of the haloes. The 
most probable Z equals to 0.7 and does not vary with ra- 
dius. The velocity ellipsoids along the halo major axis appear 
to be slightly more prolate than those in the perpendicular 
plane. This difference disappears in the innermost shell. 

4.2 Alignment 

Figure [4] shows the distribution of the angle formed between 
the major axis of the local velocity ellipsoid and the halo 
major axis (black thick profile) or the local radial direction 
(green dash-dotted profile). It is clearly seen that the local 
velocity ellipsoids in the innermost shells display a strong 
alignment with the halo major axis: the black distribution 
is peaked at 1, whereas the orientations with respect to 
the radial direction are random. The alignment is the most 
prominent in the inner part of the haloes, at r < 0.2r v , and 
gradually vanishes when approaching the virial sphere. In 
the three innermost shells, the median angle between the 
velocity ellipsoids and the halo major axis are 34, 22 and 25 
degrees. 

The alignment of the velocity ellipsoids with the halo 
major axis appears to be stronger along the major axis than 
in the perpendicular plane (see the red solid and blue dotted 
profiles in Figure fj| . The difference in degree of alignment 
between these two regions increases with radius. In both 
regions, however, the local velocity ellipsoids clearly tend to 
be aligned with the halo major axis. Exceptional deviation 
occurs around the virial sphere, where the distribution in 
the plane perpendicular to the halo major axis reveals two 
maxima associated with the alignment with the local radial 
directions and the halo major axis (see the bottom right 
panel of Figure [4} . 
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Figure 3. Triaxiality of the local velocity ellipsoids in DM haloes. 
The profiles show the distributions of the triaxiality parameter de- 
fined by the eq. J4} . Velocity ellipsoids are predominantly prolate 
at all radii. The six panels show the distributions in 6 spherical 
shells of radii as indicated in the left upper corners. The black 
(thick solid), red (solid) and blue (dotted) lines show the profiles 
for different pixel selection: all 48 pixels, 8 pixels along the halo 
major axis and 8 pixels in the plane perpendicular to the major 
axis, respectively. The black dashed profiles are the distributions 
inside relaxed haloes. 

The black dashed profiles in Figure [4] show the distri- 
butions in relaxed haloes. Compared to the sample of all 
haloes, one can see that the alignment of the local velocity 
velocity ellipsoids with the halo major axis is stronger in the 
inner parts of the relaxed haloes. On the other hand, orien- 
tations of the velocity ellipsoids at the virial sphere appear 
to be distributed in the same way. 

In terms of the alignment of the local velocity ellip- 
soids, DM haloes reveal two distinct zones: the inner part 
(r < 0.2r v ) where the ellipsoids are aligned with the halo 
major axis, and the outer part (r J> 0.5r v ) where align- 
ment with the major axis is equally probable as with the 
radial direction. The local velocity distribution in the for- 
mer is consistent with the cylindrical symmetry, whereas in 
the latter it cannot be attributed to any simple symmetry 
such as spherical or cylindrical. Interestingly, both zones co- 
incide with two characteristic regimes of the DM density 
profile: shallower and steeper than p oc r~ 2 . 

4.3 Velocity dispersion 

Figure [5] shows the distributions of the local velocity dis- 
persion in successive spherical shells. The dispersions are 
calculated using the trace of the local velocity dispersion 
tensor, i.e. a^ ot = a\ + o"| + 0"3> and are rescaled to the virial 
velocity V v = \J GM v /r v , which is a common scale of the 
velocity distributions in DM haloes of different masses. 

The distributions of the velocity dispersion have long 
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Figure 4. Alignment of the local velocity ellipsoids in DM haloes. 
The profiles show the distributions of cosines of the angle (| cosxl) 
between the major axis of the local velocity ellipsoid and the halo 
major axis (except of the green dash-dotted profiles for which 
the angle is formed between the major axis of the local velocity 
ellipsoids and the local radial direction). The distributions show 
that the ellipsoids are preferentially parallel to the halo major 
axis. Degree of the alignment increases towards the halo centres. 
The six panels show the distribution in 6 spherical shells of radii 
as indicated in the left upper corners. The black (thick solid), 
red (solid) and blue (dotted) lines show the profiles for different 
pixel selection: all 48 pixels, 8 pixels along the major axis and 8 
pixels in the plane perpendicular to the major axis, respectively. 
The black dashed profiles show the distributions inside relaxed 
haloes. 

tails at large values. These features are particularly promi- 
nent in the three outermost shells of radii r > 0.2r v . By 
comparing with the symmetric distributions representing re- 
laxed haloes (see the black dashed profiles) , we conclude that 
the high-dispersion tails are mostly populated by unrelaxed 
haloes. 

The distribution of the velocity dispersion in a given 
shell appears to be independent of the position with respect 
to the halo major axis (compare the red solid and blue dot- 
ted profiles). A small trend seems to occur in the three inner- 
most shells: the dispersion along the halo major axis tends 
to be slightly smaller than in the perpendicular plane. 

Not surprisingly, the distribution of the velocity disper- 
sion depends primarily on radius. The mean of the distri- 
bution varies with the radius from l.QV v at 0.01r v through 
1.314 at 0.1r v to 0.8V^at r v . The radial variation of the dis- 
persion distribution reflects dependance on the halo gravi- 
tational potential. 



4.4 Local anisotropy 

Here, we consider the anisotropy of the local velocity dis- 
persion tensor. The anisotropy is defined as the ratio of the 
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Figure 5. Local velocity dispersion in DM haloes, <r% t = <r? + 
o| + o| . The six panels show the distributions in 6 spherical shells 
of radii as indicated in the left upper corners. The colour coding 
is the same as in Figure [3] The distribution of the local velocity 
dispersion depends predominantly on the distance from the halo 
centres. 



Figure 6. Anisotropy 0-1/0-23 of the local velocity ellipsoids in 
DM haloes, where o| 3 = (a\ + cr|)/2. The six panels show the 
distributions in 6 spherical shells of radii as indicated in the left 
upper corners. The colour coding is the same as in Figure [4] The 
local velocity ellipsoids are highly anisotropic at all radii. Typical 
value of the anisotropy is 01/023 ~ 1.3 at radii r < 0.5r v and 
increases to 01/0-23 ~ 1.45 at the virial radius. 



major-to-minor velocity dispersions 0-1/0-23, where 0I3 = 
(erf + o|)/2. This definition ignores the fact of triaxiality 
of the velocity ellipsoids and, therefore, provides only a par- 
tial description of the ellipsoids. Nevertheless, keeping in 
mind that the local velocity ellipsoids are preferentially pro- 
late, i.e. (Ti — 02 > 02 — 03, we expect that the 0-1/0-23 
ratio provides the first order and single-parameter descrip- 
tion of the local velocity distribution (we will consider both 
the 01/02 and 01/0-3 ratios in the following section, where 
we compare the local and spherically averaged properties of 
the velocity dispersion tensor). Contrary to the definition of 
the anisotropy parameter /3, here we assume the cylindrical 
symmetry of the local velocity distribution. 

Figure [6] shows the distribution of the 01/0-23 ratio in 6 
spherical shells. The distributions inside the innermost shells 
of radii r < 0.2r v depends barely on radius. The mode is 
equal to 1.3 in the first shell and 1.25 in the following three 
shells. The distributions have long tails extending to high- 
value anisotropies. Similar to the case of the velocity disper- 
sion, the unrelaxed haloes appear to populate preferentially 
the tails of the distributions (compare the black thick and 
dashed profiles corresponding to the all and relaxed haloes). 
The anisotropy along the major axis of the virial sphere 
tends to be smaller than in the perpendicular plane. 

The local velocity ellipsoids at the virial sphere tend to 
be substantially more anisotropic than in the inner parts. 
The mode occurs at 01/023 = 1.45 and the distribution has 
a slowly decreasing tail at high values. Distributions for re- 
laxed and unrelaxed haloes appear to be undistinguishable. 
The anisotropy around the halo major axis in these shells 
tend to be smaller than in the perpendicular plane, in con- 



trast with the innermost shells (compare the red solid and 
blue dotted profiles). 

Similar to the alignment of the velocity ellipsoids, sta- 
tistical properties of the 01/023 ratio seem to differentiate 
the central parts of the haloes from those around the virial 
sphere. The former are characterised by self-similar distri- 
butions with the maximum at approximately 1.25, whereas 
the latter by the distribution gradually getting wider and 
shifted to higher values when approaching the virial sphere. 
The transition between these two zones is continuous. The 
approximated radius of the transition is 0.2r v at which the 
shape of the distribution starts to deviate from that found 
in the shells of small radii (see the middle right panel of 
Figure . 



5 SPHERICALLY AVERAGED ANISOTROPY 

In this section, we address the problem how to measure 
the anisotropy of the velocity dispersion tensor in spherical 
shells. It is naturally expected that the velocity dispersion 
computed in a shell should represent information contained 
in the distribution of the local properties. In particular, 
spherically averaged profiles (computed in spherical shells) 
are expected to recover the most probable values of the local 
counterparts. In order to achieve this, the velocity dispersion 
tensor should conform with with the symmetry preferred by 
the local velocity ellipsoids. Strong alignment of the velocity 
ellipsoids with the halo major axis implies that the velocity 
dispersion tensor should be calculated using a Cartesian co- 
ordinate system or a cylindrical one aligned with the halo 
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Figure 7. Radial profiles of the anisotropy of the velocity disper- 
sion tensor calculated in spherical shells. The red solid and blue 
dotted lines are the profiles for a\/o2 and (71/0-3, respectively 
(with dashed lines showing the first and third quartiles). The 
points with error bars represent the distributions of the local val- 
ues measured in 7 spherical shells (the median and the error bar 
corresponding to the 50 per cent probability range). The <t\/<T2 
and <ti /<T3 ratios computed in spherical shells reproduce the local 
values. Discrepancy at radii r > 0.5 V occurs due to more random 
orientations of the velocity ellipsoids in this part of the haloes. 

axis (if one ignores triaxiality of the ellipsoids). Here we 
follow the former approach and calculate the velocity dis- 
persion in spherical shells using the eq. ([3|. 

Figure shows the radial profiles of the anisotropy of 
the velocity dispersion tensor calculated in spherical shells 
around the centres of all haloes. We plot both the ratio of the 
major-to-median (01/02) and the major-to-minor (0-1/03) 
velocity dispersions. We used the same shells in all haloes. 
The profiles are the lines connecting the median values mea- 
sured in every shell and the errors are given by the first and 
third quartiles. 

The spherically averaged profiles are compared with the 
the distribution of the local values measured in the shells 
defined in the previous section. The distributions are rep- 
resented by the quartiles (black points with the error bars 
showing the median and the 50 per cent probability range). 

The spherically averaged profiles of the anisotropy are 
fairly consistent with the local values at radii r < 0.5r v . 
The quartile profiles of the major-to-median velocity dis- 
persion ratio, 01 /(T2, trace the quartiles of the local values, 
whereas the 0-1/0-3 profiles tend to slightly underestimate 
the local counterparts. In this range of radii, the anisotropy 
is well-approximated by a flat profile with 01/0-2 ~ 1-25 and 
0-1/0-3 » 1.42. 

The spherically averaged anisotropy appears to under- 
estimate the true local values at radii r > 0.5r v . This is par- 
ticularly well visible for the 0-1/03 ratio, for which deviation 
between the median of the local and spherically averaged 
values of the anisotropy is of the order of 2o at r ~ r v . 
This discrepancy is mostly caused by a tension between the 
assumed symmetry of the velocity ellipsoids and the true 
orientations of the local velocity ellipsoids. Larger scatter in 
the alignment of the local velocity dispersions (more ran- 
dom orientations) produces artificially more isotropic veloc- 
ity dispersion tensor calculated in spherical shells. 

Figure [8] shows the profiles of the angle between the 
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Figure 8. Alignment of the velocity ellipsoid calculated in spher- 
ical shells. The lines show the 50 (upper) and 68 (lower) per cent 
probability of the angle (cos \) between the major axis of the ve- 
locity ellipsoid measured in spherical shells and the halo major 
axis. The velocity ellipsoids representing the velocity dispersion 
tensor calculated in spherical shells are aligned with the halo ma- 
jor axis. The red solid and black dashed lines show the profiles 
for DM particles and subhaloes, respectively. 



major axis of the velocity ellipsoid measured in spherical 
shells and the halo major axis. The lines indicate the 50 and 
68 per cent probability range of the angle measured at every 
radius (with the most probable values at all radii equal to 
|cosx| = !)■ This clearly shows that the velocity ellipsoids 
computed in spherical shells are aligned with the halo major 
axis. The scatter of the angles is comparable to the scatter 
of the local values (compare with Figure [3} . 

For the sake of illustration how /3 parameter misrepre- 
sents the true picture of velocity anisotropy in DM haloes, 
we compare in Figure [9] the spherically averaged profiles of 
P and its analogue measured in a cylindrical coordinates 
system, i.e. 1 — (erf + °'I)/(2o'i)- As expected, the j3 pro- 
file differs substantially from its counterpart calculated in 
cylindrical coordinates. In particular, nearly isotropic veloc- 
ity distribution (in spherical coordinate system) in the halo 
centre, i.e. /3 « 0, is an effect of angular averaging over veloc- 
ity ellipsoids which are strongly aligned with the halo major 
axis. The true orientations of the local velocity ellipsoids vi- 
olate some crucial assumptions underlying the definition of 
the anisotropy parameter j3 (radial orientations of the lo- 
cal velocity ellipsoids). Violation of these assumptions leads 
to artificially more isotropic velocity dispersion tensor (in 
terms of /3). 

In Figure |§J we also compare the spherically averaged 
profiles with the local values. The anisotropy defined in 
cylindrical symmetry recovers the local values at r < 0.5r v . 
As in the case of Figure [7J the discrepancy at large radii is 
due to more random orientations of the local velocity ellip- 
soids. On the other hand, the anisotropy parameter /3 gives 
erroneous impression of significantly more isotropic veloc- 
ity distribution. The median /3 profile lies (2 — 4)o lower 
than the median profile of the true local anisotropy. Using 
j3 parameter leads to artificial isotropisation of the velocity 
distribution not only in the centres of DM haloes, but also 
at radii comparable to the virial radius. 

The anisotropy parameter j3 leads to a misleading pic- 
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Figure 9. Radial profile of the anisotropy of the velocity dis- 
persion tensor calculated in spherical shells. The blue dotted and 
red solid profiles show the anisotropy parameter /3 = 1 — (<t| + 
t7^)/(2cr^) and its counterpart defined in cylindrical coordinate 
system, i.e. 1 — (<r| + ct|)/ (2<rf). Black points with error bars rep- 
resent the distributions of the local values measured in 7 spher- 
ical shells (the median and the error corresponding to the 50 
per cent probability range). The anisotropy defined in cylindri- 
cal coordinate system and computed in spherical shells recovers 
the local values. Small discrepancy at radii r > 0.5 V is due to 
more random orientations of the velocity ellipsoids in this part of 
the haloes. The anisotropy parameter is inconsistent with the 
true anisotropics of the local velocity ellipsoids. Its significantly 
smaller values result from averaging axially symmetric velocity 
ellipsoids in spherical coordinate system. The black dash-dotted 
line shows the median 1 — <r| 3 /cj profile calculated with the use 
of subhaloes. Dashed lines show the first and third quartilcs. 



ture of the true anisotropy of the velocity dispersion ten- 
sor not only when applying to spherical shells, but also lo- 
cally. For example, the median f3 in the equatorial pixels of 
the innermost shells (see Figure [T]) point to nearly isotropic 
or weakly tangentially-biased velocity dispersion tensor, i.e. 
/3 <J 0. In fact, this is a consequence of a peculiar orientation 
of the local velocity ellipsoids whose major axes are prefer- 
entially parallel to the direction of the polar angle in the 
equatorial stripe. Averaging over the angles, a commonly 
used routine based on adopting the mean (cr| + a^)/2 as 
the variance of the tangential velocity, mixes the major and 
minor axes of the velocity ellipsoids and thus gives rise to 
an artificial isotropisation. 



6 DISCUSSION AND CONCLUSIONS 

We studied the properties of the local velocity dispersion 
tensors in cluster-size simulated DM haloes. We found that 
the velocity ellipsoids representing the tensors are strongly 
aligned with the halo major axis defined as the axis minimis- 
ing the moment of inertia calculated inside the virial sphere. 
Statistical properties of the orientations and anisotropies of 
the local velocity dispersion tensor do not vary with radius 
in the central parts of the haloes (r < 0.2r v ). At large radii, 
the orientations become gradually randomised and the ellip- 
soids more elongated. These two distinct zone of the haloes 
coincide with two characteristic regimes of the DM density 
profile: shallower and steeper than p oc r~ 2 . 

The following simplified model captures the main fea- 
tures of the velocity ellipsoids. In the reference frame whose 
axes are aligned along the principle axes of the density dis- 
tribution (with x, y and z axes corresponding to the major, 
medium and minor axis, respectively), the number-density n 
as the function of velocities V = {v x ,v y ,v z } and coordinates 
r = {a;, y, z} can be written as : 



n(V,f) 

P y (x,y,z) 
Pz{x,y,z) 



no(r) 



37T 3 cr|crgcrf ) x / 2 



exp 



2cr| 2cr2 2a1 



wl(0A + w) 



1 - -| = p 0z 1 



V 



1+w 
1+w 



+ 



y_ 

"sy 



+ 



(5) 

(6) 
(7) 
(8) 

(9) 



where no(x, y, z) is the density profile of the halo and r ax ,r ay 
and r sz are the sca le radii of the NFW density profile 
l)Navarro et al.l Il997l ) fitted along the major, medium and 
minor axis, respectively. For cluster-size haloes studied in 
the paper 



0.35, p 0z 



0.45 



(10) 



5.1 Subhaloes 

The number of subhaloes is too small to study the proper- 
ties of the local velocity dispersion tensor. Here, we calcu- 
late spherically averaged profiles of the anisotropy and the 
alignment, and compare with those obtained for DM parti- 
cles. We use subhaloes found with the BDM algorithm and 
containing at least 30 DM particles. 

The black dashed lines in Figure [8] and Figure [9] show 
the profiles calculated with the subhaloes (assuming equal 
weights). In the range of radii at which subhaloes are de- 
tectable, the anisotropy of the velocity ellipsoids does not 
differ from that for DM particles. The velocity ellipsoids are 
also aligned with the halo major axis, although the scatter 
appears to be slightly larger (see Figure [8]). 



Alignment of the local velocity ellipsoids is inconsis- 
tent with spherically symmetric velocity distribution (radi- 
ally oriented ellipsoids) that is an assumption underlying 
definition of the anisotropy parameter /3. As a consequence, 
using the anisotropy parameter j3 (assuming radial orienta- 
tions of the local velocity ellipsoids) leads to an erroneous 
picture of significantly more isotropic velocity distributions 
than they really are. Typical ratio of the major-to-minor 
axis of the local velocity ellipsoids is equal to 1.3 at radii 
r < 0.5r v and increases strongly with radius in the outer 
part of the haloes. The spherically averaged (calculated in 
spherical shells) profiles of the anisotropy reproduce the lo- 
cal values when one uses a Cartesian or cylindrical coordi- 
nate system that respects the symmetry of the alignment of 
the velocity ellipsoids with the halo major axis. 
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6.1 Phase-space properties of DM haloes 

The orbital anisotropy is an important diagnostic in study- 
ing the phase-space properties of DM haloes. Contrary 
to the previous find ings based on the a nisotropy pa- 
rame t er P (see e.g. Diemand et al.l 20041; Woitak et alj 



20051; lAscasibar fc Gottlobej|2008l ; llannuzzi fc Dolad 12012} 
Lemze et all 120121 1 . our study shows that the velocity dis- 
tribution in the most relaxed parts of the haloes is highly 
anisotropic and the local velocity ellipsoids are aligned with 
the major axis of the halo shape (cylindrical symmetry). 
Bearing in mind that DM haloes are aligned with the most 
prominent filaments determining pr eferential directions o f 
major mergers and smooth accretion (|Libeskind et al J2012t ). 
this implies that the orbital structure in DM haloes is likely 
related to the formation processes of the haloes. 

Further studies on a connection between the phase- 
space structure of cosmological haloes and their formation 
history should drop the assumption of spherical symmetry. 
This is also important in addressing such problems as to 
what extent and in what sense cosmological haloes are re- 
laxed objects. In the light of the results presented in this 
paper, the assumption of spherically symmetric velocity dis- 
tributions in DM haloes acts as an artificial phase mixing 
which may lead to a picture, in which DM haloes appear to 
be more equilibrated than they are. 



6.2 Dynamical modelling of quasi-spherical 
systems 

Although our studies are based on the specific sample of 
DM haloes (cluster-size haloes), the overall picture of the 
velocity ellipsoids aligned with the halo major axis may 
likely be generic. The argument supporting this expecta- 
tion is twofold. First, the alignment is more prominent in 
relaxed haloes rather than in recent mergers. Second, all 
processes of anisotropic halo formation are common to all 
DM haloes and, therefore, the same mechanisms breaking 
spherical symmetry of the velocity distribution should oper- 
ate at all scales. 

Our results may have some consequences for the mass 
modelling of kinematical data in such objects as galaxy clus- 
ters, elliptical galaxies and dwarf spheroidals. Most of dy- 
namical models assume spherically symmetric velocity dis- 
tributions whose anisotropy is quantified with the /? parame- 
ter. It has barely ever been verified whether this assumption 
conforms with the true orbital s tructure inside th e se obj ects. 
This problem was addressed bv lCappellari et al.l (|2007T ) in a 
detailed study of the stellar kinematics in massive ellipticals. 
It was shown that the velocity ellipsoids are preferentially 
aligned with the galaxy shape, on the contrary to what is 
assumed when considering /3 parameter. Some signatures of 
similar orientation of the veloc ity ellipsoids in galax y clus- 
ters were also demonstrated bv lSkielboe et all |2012T ). 

Alignment of the local velocity ellipsoids with the ma- 
jor axis may inevitably affect the mass inference with the 
use of dynamical models assuming radial orientations of the 
ellipsoids. An additional analysis is required to address this 
question in a fully quantitative way. However, in order to in- 
troduce the scale of the problem, we note that, in the light 
our results, the projected velocity dispersions in two perpen- 
dicular directions may differ by as much as 50 per cent (see 



Figure [7]). Using the M v oc a 3 scaling relation, this leads to 
the ratio of the mass estimates in both projections equal to 
3.4. 



6.3 Observational constraints on /3 

The anisotropy parameter /3 has been measured for a 
number of astrophysical systems. The estimates obtained 
from the stacked k i nema t ical data of galaxies i n clusters 
(IBiviano fc KatgertJ I200I ; IWoitak fc Lokasl 120101) and the 
satell ite galaxies around isolated hosts ( Woitak fc Mamonl 



are fairly compatible with the spherically averaged 
profiles from the simulations. We emphasise, however, that 
this consistency does not imply the fact that /3 describes the 
true orbital anisotropy, because stacking kinematical data 
introduces artificial sphericalisation of the phase-space dis- 
tribution that is equivalent to measuring /3 in spherical shells 
of simulated objects. 

Observational constraints on the spherical anisotropy 
j3 in individual clusters exhibit a huge variety of profiles 
llBenatov et~aH 120061; IHwang fc Lee1l200g| ; IWoitak fc Lokasl 
l20ld ; IWoitak et al.ll2009F ~ Some measurements are far out- 
side the margins allowed by t he spherically averag ed profiles 
of simulated clusters (see e.g. IHwang fc Leg|2008l ). We sus- 
pect that this substantial scatter of observational constraints 
on /3 in individual clusters may result from ignoring the angle 
between the line of sight and the major axis of the veloc- 
ity ellipsoids. This projection effect needs to be addressed 
in future in a quantitive analysis of mock kinematical data 
generated from cosmological simulations. 
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